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Cerebral blood flow identifies responders to 
transcranial magnetic stimulation in auditory verbal 
hallucinations 

P Homan 1 ' 3 , J Kindler 1 ' 3 , M Hauf 2 , D Hubl 1 and T Dierks 1 

Auditory hallucinations comprise a critical domain of psychopathology in schizophrenia. Repetitive transcranial magnetic 
stimulation (TMS) has shown promise as an intervention with both positive and negative reports. The aim of this study was to test 
resting-brain perfusion before treatment as a possible biological marker of response to repetitive TMS. Twenty-four medicated 
patients underwent resting-brain perfusion magnetic resonance imaging with arterial spin labeling (ASL) before 10 days of 
repetitive TMS treatment. Response was defined as a reduction in the hallucination change scale of at least 50%. Responders 
(n = 9) were robustly differentiated from nonresponders (n = 15) to repetitive TMS by the higher regional cerebral blood 
flow (CBF) in the left superior temporal gyrus (STG) (P<0.05, corrected) before treatment. Resting-brain perfusion in the left 
STG predicted the response to repetitive TMS in this study sample, suggesting this parameter as a possible bio-marker of 
response in patients with schizophrenia and auditory hallucinations. Being noninvasive and relatively easy to use, resting 
perfusion measurement before treatment might be a clinically relevant way to identify possible responders and nonresponders 
to repetitive TMS. 
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Introduction 

Hallucinations are perceptual experiences in the absence of 
an adequate external stimulus, and auditory verbal hallucina- 
tions (AVH) are a key symptom of psychosis. In schizo- 
phrenia, the 1 -month prevalence of these hallucinations 
exceeds 70%, 1 and in 25-30% of patients these perceptions 
are resistant to medication, resulting in functional disability 
and a low quality of life. 2,3 With limited therapeutic options, the 
development of new treatment strategies 4 is needed. One 
such intervention is repetitive transcranial magnetic stimula- 
tion (TMS). Hoffman was the first to report that low-frequency 
TMS decreased the severity of medication-resistant AVH 
when the coil is directed at the left temporoparietal 5 area. 
Since that initial report, a number of studies have been 
conducted with a similar treatment paradigm and focus. 6 " 21 
Nevertheless, in a recent study with the largest sample size so 
far, no difference between sham and stimulation has been 
found 22 However, it has been suggested that patients 
respond differentially to TMS, 23 making the exploration for 
possible predictors desirable. Predictors may be either found 
in the clinical phenomenology of symptoms or in brain 
function. With respect to symptoms, only minimal response 
to TMS has been reported in continuous hallucinators. 6 ' 14 
With regard to brain function, only functional magnetic 
imaging (fMRI), but not resting perfusion, has been used to 
predict response to TMS in AVH. Using temporoparietal TMS, 
lower levels of hallucination-related activation in Broca's area 



predicted the response rate in intermittent hallucinators, and 
reduced coupling between Wernicke's and a right homolog of 
Broca's area predicted the response-rate in continuous 
hallucinators. 15 fMRI requires detecting the symptoms of 
hallucinations in the scanner, which is extremely challenging 
depending on many uncontrollable factors, thus making it 
clinically practical only in a few patients. Therefore, predictors, 
such as resting neuronal activity, not requiring symptom and 
patient cooperation, would clearly be clinically more useful. 
We recently showed that the improvement of symptoms after 
10 days of TMS treatment was associated with decreased 
neuronal activity in the primary auditory cortex, Broca's area 
and cingulate gyrus 24 suggesting that resting perfusion, 
reflecting neuronal activity, is an appropriate marker for the 
effectiveness of TMS. In the current study, we applied a 
different approach. Our aim was to identify responders and 
nonresponders post focusing resting-state perfusion and not 
to test the effectiveness of TMS itself. Considering that there 
is a hyperactivity in language-related regions in AVH patients, 
which is related to symptom severity, 25-27 we postulated that 
cerebral-resting perfusion before TMS treatment is higher in 
responders compared with nonresponders. To achieve a 
larger sample size, we here pooled patients using two 
different TMS protocols. A repetitive TMS protocol that has 
recently been introduced into clinical research is theta-burst 
stimulation. 28-30 Compared with the often-used 1-Hz stimula- 
tion protocol, theta-burst stimulation has the advantage that 
the application duration is very short, making the treatment 
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more convenient for the patient. In a recent study, we found 
theta-burst and 1-Hz TMS to be equally effective, 24 and 
therefore we randomly assigned patients to receive either a 1 - 
Hz or theta-burst TMS protocol. Before treatment, we 
measured whole-brain resting perfusion using magnetic 
resonance (MR) arterial spin labeling (ASL), an MR technique 
that provides a direct quantitative measure of cerebral blood 
flow (CBF). 31-33 ASL is a noninvasive technique that reliably 
measures whole-brain CBF with converging results to those 
obtained with invasive PET(Positron emission tomography) 
perfusion imaging. 34 

Materials and methods 

Patients and clinical investigation. Inclusion criteria were 
the diagnosis of schizophrenia or schizoaffective disorder 
according to ICD-10 (F20, F25), age between 18 and 65 
years, right-handedness and medication-resistant AVH. 
Therapy refractoriness was defined as having no response 
to at least two different antipsychotic treatment options within 
licensed dose ranges, each administered for at least 8 
weeks. Medication dose remained stable during the treat- 
ment. Exclusion criteria were history of epileptic seizures, 
signs of elevated neuronal activity by electroencephalogra- 
phy, MR contraindications and medical disorders other than 
schizophrenia or schizoaffective disorder. Substance misuse 
in the 4 weeks before treatment was ruled out by a urinary 
drug screen before treatment. All patients underwent 
identical diagnostic procedures and identical MR protocols. 
The diagnostic procedure was performed on the basis of 
clinical interviews and psychiatric history. The assessment of 
psychopathology, consisting of the auditory hallucination 
rating scale (AHRS), 35 the hallucination change scale 35 and 
the positive and negative symptom scale (PANSS), 36 was 
performed before and after 10 days of TMS treatment. To 
detect clinically relevant improvements in symptoms, we 
used a conservative level of at least 50% in the hallucination 
change scale to define response to TMS, similar to the one 
used in previous studies. 14 ' 35 The Edinburgh handedness 
scale 37 was assessed to ensure that all patients were right- 
handed. The investigation was conducted in accordance with 
the Declaration of Helsinki and approved by the local ethics 
committee. Patients provided informed written consent to 
participate in the study. 

Experimental design. MRI was conducted on a 3-T whole- 
body MRI system (Magnetom Trio, Siemens Medical 
Systems, Erlangen, Germany) with a standard 12-channel 
radiofrequency head coil. High-resolution three-dimensional 
structural (s)MRI and ASL were acquired before TMS. T1- 
weighted three-dimensional magnetization-prepared rapid 
gradient echo scans were recorded (number of slices = 176, 
matrix = 256 x 256, slice thickness = 1 mm, voxel size = 1 
x 1 x 1 mm 3 ), serving as a high-resolution three-dimensional 
anatomical template for coregistration with functional data. A 
pseudo-continuous ASL technique was used to measure 
CBF. 38 In this gradient-echo echo-planar imaging sequence, 
interleaved images with and without labeling were acquired. 
A delay of 1250 ms was applied between the end of the 
labeling pulse (label time = 1720 ms) and image acquisition 



(slice acquisition time = 45 ms), in order to reduce transit 
artifact (field of view = 220 mm 2 , matrix = 64 x 64, repetition 
time/echo time = 4000/1 8 ms and flip angle = 25°). A total of 
14 slices (voxel size = 3.4 x 3.4 x 6 mm 3 , gaP=1.5mm) 
were acquired in the anterior and posterior commissure line, 
from inferior to superior, in a sequential order. The pseudo- 
continuous ASL scan comprised 80 acquisitions. ASL data 
analysis was performed using aslm 39 (downloadable at 
http://aslm.sourceforge.net), a newly developed toolbox that 
allows multiple subjects to be analyzed at once and to 
perform region of interest analysis of ASL data, aslm is an 
object-oriented toolbox based on Matlab (Matlab version 8, 
release 14; The MathWorks, Natick, MA, USA), SPM 8 
(Wellcome Department of Imaging Neuroscience, London, 
England; http://www.fil.ion.ucl.ac.uk/spm/software/spm8/) and 
marsbar. 40 All ASL time series were first realigned to correct 
for motion artifacts. We calculated flow-time series by subtract- 
ing the labeling images from the control images, and sub- 
sequently computed mean CBF images for each subject. 41 
Each individual subject's T1 anatomy was segmented into 
gray matter (GM) and white matter. Mean ASL images were 
then coregistered to GM-segmented T1. T1, GM, white matter 
and ASL images were normalized to the SPM MNI T1 template. 
ASL images were spatially smoothed with a three-dimensional 
8-mm Full Width at Half Maximum Gaussian kernel. Data were 
z-transformed (z= (voxel CBF-global GM CBF)/s.d.) to remove 
sources of variance caused by differences of global mean CBF 
between acquisitions and corrected for GM using GM 
segments as inclusive masks. For visualization, CBF maps 
and clusters of the whole-brain analysis were overlaid onto 
T1 -weighted MR image slices using aslm. 

Neuronavigated TMS. We chose area Spt (Sylvian parie- 
totemporal), located in the sylvian fissure at the parietotem- 
poral boundary, as the target region for TMS. Area Spt is a 
sensorimotor interface between the sensory and motor 
speech systems, 42 is strongly left dominant, 43 and is 
activated equally by the perception and reproduction of 
aurally or visually presented words. There is evidence of 
functional abnormalities within area Spt in patients with 
schizophrenia. 44 Because area Spt is connected to the 
frontal and temporal language regions that are involved in the 
generation of AVH, it qualifies as an ideal target region for 
TMS modulation of AVH. A modified version of an evaluated 
language-processing task introduced by Hickok et al. 43 was 
used to localize the target region (Spt) for TMS treatment in 
an fMRI scan of each patient in the first MR measurement. 
This measurement took place immediately after the resting 
perfusion measurement with ASL. After superimposing the 
individual's fMRI activation map on its anatomical mesh 
(three dimensional reconstruction of T1 image), a TMS target 
point was created at area Spt. A custom TMS stimulator 
(MagPro R 100, Medtronic Functional Diagnostics, Skov- 
lunde, Denmark) was used to generate repetitive biphasic 
magnetic pulses. Magnetic pulses were delivered with a 
figure-8 coil (Magnetic Coil Transducer MC-B70, Medtronic). 
Stimulation of the motor cortex with single TMS pulses until a 
movement of the contralateral thumb was detected permitted 
the identification of the individual resting motor thresholds. 45 
The center of the coil was held tangentially to the skull. TMS 



Translational Psychiatry 



Response to transcranial magnetic stimulation 

P Homan et al 



Table 1 Patient characteristics with mean values, s.d., statistical differences with P-values 



Characteristic 


Responders (n = 9) 


Non-responders (n = 15) 


Test 


Value 


df 


P 


Gender, no. F/M 


4/5 


10/5 


Fisher's exact 


NA 


NA 


0.4 


Diagnosis, no. F20/F25 


7/2 


14/1 


Fisher's exact 


NA 


NA 


0.5 


Age, mean (s.d.), years 


38.3(13.9) 


44.4(10.1) 


f-test 


1.3 


22 


0.2 


Age at onset, mean (s.d.), years 


24.7 (6.5) 


26.6 (6.6) 


Mest 


-1.4 


22 


0.2 


Duration of disease, mean (s.d.), years 


12.9 (7.1) 


18.6(10.6) 


f-test 


-1.2 


22 


0.2 


Chlorpromazine equivalent dose at study entry, mean (s.d.) 


448.9 (132.2) 


586.2 (257.1) 


f-test 


-1.5 


22 


0.1 


PANSS score before TMS, mean (s.d.) 


75.8 (17.7) 


72.6 (15.4) 


f-test 


0.5 


22 


0.6 


AHRS score before TMS, mean (s.d.) 


33 (6.7) 


35.7 (6.1) 


f-test 


1.0 


22 


0.3 


Global mean CBF, corrected for GM(ml 100g" 1 min" 1 ) 


58.3(12.5) 


52.0 (9.0) 


f-test 


1.5 


22 


0.2 



Abbreviations: AHRS, auditory hallucination rating scale; CBF, cerebral blood flow; df, degrees of freedom; F, female; F20, schizophrenia; F25, schizoaffective 
disorder (according to ICD-1 0); GM, gray matter; M, male; NA, not applicable; PANSS, positive and negative syndrom scale; TMS, transcranial magnetic stimulation. 



pulse intensity was adjusted to 90% of the resting motor 
threshold. Patients were randomly assigned to receive a 1- 
Hz (n=12) or theta-burst (n=12) TMS protocol. In both 
treatment groups, the target area was stimulated for 10 
consecutive days. Stimulation at 1 Hz was applied once a 
day, according to the protocol used by Hoffman etal. 35 In the 
theta-burst group, each burst contained three pulses at 
30 Hz, repeated with an interburst interval of 100 ms. 29 Theta 
burst was applied in double trains with a 15-min intertrain 
interval. On the first 3 days, 2 double trains of theta burst 
were applied (based on the protocol in Nyffeler et al. 29 ), 
whereas on days 3-10 one double-theta-burst train was 
applied. Safety protocols were in accordance with interna- 
tional safety standards of TMS experimentation. 46 A frame- 
less, ultrasound-based stereotactic system was used for 
neuronavigation (BrainvoyagerTMS Neuronavigator System, 
Brain Innovation, BV 2006, 47 Maastricht, The Netherlands). 
After superimposing the individual's fMRI activation map on 
its anatomical mesh (three-dimensional reconstruction of T1 
image), a TMS target point was created at area Spt. Details 
on the fMRI stimulation paradigm have recently been 
published elsewhere. 24 

Statistical analysis: CBF. Average mean GM CBF values 
of responders and nonresponders were compared in a two- 
sample Mest. A global and assumption-free second-level 
random-effect analysis was then performed using a two- 
sample Mest on GM CBF images obtained in each 
participant (degree of freedom (df) = 22). We report clusters 
with a voxel-level threshold of P< 0.001 (uncorrected) at a 
cluster-level threshold of P<0.05, whole-brain corrected for 
familywise error. 

Results 

Clinical data. Twenty-four patients were included in the 
study. Clinical and demographic characteristics of the subject 
sample are detailed in Table 1. After 10 days of TMS 
treatment, patients were rated again by a blinded clinician. 
Nine patients showed a clinical improvement of at least 50% 
and were therefore classified as responders; the other 15 
patients were classified as nonresponders. Responders and 
nonresponders did not differ in the levels of AHRS and 
PANSS before treatment (Table 1). In addition, the frequency 
item (P=0.8), the reality item (P=0.4) and the loudness 
item (P=0.1) of the AHRS did not differ in responders and 



nonresponders. Responders and nonresponders were not 
different in age, gender, age of onset, disease duration or 
medication levels (Table 1). No differences were detected 
between the 1-Hz (mean: 0.32 ±0.39) and theta-burst 
(mean: 0.32 ±0.32) protocols comparing the amount of 
improvement as measured with the hallucination change 
scale (P= 1). There was no difference in the relationships of 
responders to nonresponders in the 1-Hz group (4:8) and in 
the theta-burst group (5:7, P= 1). 

ASL data. The average mean global CBF in the responders 
and nonresponders was 58.3 (±12.5) and 52.0 
(±9.0)ml100g~ 1 min -1 , respectively (Table 1). The aver- 
age mean global CBF images of responders and nonrespon- 
ders are shown in Figure 1 . As average mean global CBF did 
not differ between responders and nonresponders (P=0.2), 
analyses were performed with normalized data to reduce 
variability introduced by nonspecific fluctuations in global 
CBF. In the whole-brain random-effect analysis, two contrasts 
were tested: first, responders > nonresponders, that is, voxels 
that were significantly higher perfused in responders than in 
nonresponders before treatment, and second the inverse 
contrast (nonresponders > responders). In the first contrast, 
one single cluster of 120 voxels survived a cluster-level 
whole-brain correction for familywise error in the left superior 
temporal gyrus (STG) (Brodmann area 22, peak: t [22] = 7.28, 
x/y/z= -58/- 26/4, Figure 2). No clusters survived in the 
inverse contrast (nonresponders > responders). 

ASL data and psychopathology. The normalized levels of 
CBF in the left STG before treatment were correlated 
(Spearman's rank correlation) with the amount of improve- 
ment, as measured with the hallucination change scale 
(r=0.7, df: 22, P< 0.0001, Figure 3). There was no relation- 
ship between the level of perfusion in STG in responders and 
nonresponders and the levels of symptoms as measured by 
the PANSS (r= -0.04, df: 22, P=0.9) before treatment. In 
addition, pretreatment scores of the AHRS items frequency 
(P=0.9), reality (P=0.7) and loudness (P=0.5) were not 
correlated with the level of STG perfusion. 

Reciever operator characteristic (ROC) analysis. To 

further investigate the discriminative power of CBF in the 
left STG, we conducted a receiver operator characteristic 
(ROC) analysis. That way, we measured how well CBF in the 
left STG could be used to distinguish between responders 
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Figure 1 Average mean global CBF images of responders (top) and 
nonresponders (bottom). CBF values were z-transformed. 
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Figure 2 Contrast responders > nonresponders, that is,voxels that were 
significantly higher perfused in responders than in nonresponders before treatment. 
Single cluster of 120 voxels survived a cluster-level whole-brain correction for 
familywise error in the left superior temporal gyrus (Brodmann area 22, peak: 
t[22] = 7.28, x/y/z = - 58/ - 26/4). No clusters survived in the inverse contrast 
(nonresponders > responders). 



and nonresponders. The area under the ROC curve (AUC) 
was 0.96 (s.e.: 0.05, 95% confidence interval: 0.87-1.0, 
P< 0.0001 , Figure 4). The cutoff point for best sensitivity and 
specificity was calculated at z=0.26 (Figure 4). 

Discussion 

Here we could demonstrate that CBF predicts TMS ther- 
apeutic response in patients with schizophrenia suffering from 
AVH. Responders showed a higher resting perfusion in the left 
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Figure 3 Scatter plot of improvement in percent, as measured with the 
hallucination change scale and CBF in the left superior temporal gyrus before 
treatment (Spearman's rank correlation, r= 0.7, df: 22, P< 0.0001). CBF values 
were z-transformed. 
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Figure 4 ROC analysis to test the discriminative power of the CBF in the left 
superior temporal gyrus before treatment in responders and nonresponders to 
transcranial magnetic stimulation. The AUC was 0.96 (s.e.: 0.05, 95% confidence 
interval: 0.87-1.0, P< 0.0001). The cutoff point for best sensitivity and specificity 
was calculated at z=0.26. 



STG before treatment. The clinical relevance of this finding 
derives from the fact that possible responders to TMS might 
be identified before the treatment by a single noninvasive MR 
measurement of resting perfusion. 

In the present study, a sample of 24 patients with AVH who 
had been treated for 10 days with TMS was examined to test 
whether resting perfusion measurement with ASL before 
treatment could be used to distinguish between responders 
and nonresponders. A single cluster of 120 voxels in the left 
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STG was detected that was significantly higher perfused in 
responders than in nonresponders and that was very close to 
the region targeted by TMS, suggesting that the efficacy of 
TMS depends on the level of brain activation in the area close 
to the target region. Patients with higher brain activity in the 
left STG might constitute a specific clinical subgroup of 
patients responsive to TMS. This is in line with recent findings. 
The reduced but still significant effect size of TMS on 
refractory AVH, when compared with previous quantitative 
reviews that has recently been published, 23 suggests the 
possibility that TMS may be effective in clinical subgroups of 
patients only. In addition, a study with the largest sample size 
(n = 62) so far showed that low-frequency TMS administered 
to the left temporoparietal cortex or to the site of maximal 
hallucinatory neuronal activation is not more effective for 
medication-resistant AVH than sham treatment. 22 Although 
the authors reported that mean severity of AVH significantly 
decreased over time, no significant difference in reduction of 
AVH or any other psychotic symptom was revealed among the 
three treatment conditions. They described that even when 
the two groups receiving guided and nonguided TMS were 
combined (n = 42), TMS was not superior to sham treatment 
in reducing AVH or other psychotic symptoms. However, 
resting brain perfusion was not measured in their study, and 
thus the missing effect might be explained by lower levels of 
CBF corresponding to a low neuronal activity in the target 
region. On the other hand, in the aforementioned study, an 
online fMRI paradigm was used during the experience of AVH 
to determine the TMS target region. As BOLD signal changes 
represent relative changes of neuronal activity, they may be 
associated with unresponsiveness to TMS. Thus, the present 
study extends the previous ones by suggesting that baseline 
resting CBF in specific language-related cerebral areas is the 
key factor determining the efficacy of TMS therapy. More 
specifically, a nonresponder showing statistically relevant 
BOLD activations during AVH might still be 'under-perfused' in 
the same area compared with a responder. Our finding of an 
elevation in the left STG perfusion in responders to TMS gives 
further evidence for the involvement of primary language 
production and perception areas in AVH. It is widely believed 
that AVH arise from a disorder of inner verbal experiences and 
its monitoring. 48 ' 49 Altered activation in speech production 
areas (that is inferior frontal gyrus), and altered coupling with 
monitoring areas (anterior cingulate) and language reception 
areas (Wernicke's area), has been proposed in a model of 
bottom-up dysfunction through overactivation in the second- 
ary (and occasionally primary) sensory cortices, which leads 
to the experience of vivid perceptions in the absence of 
sensory stimuli. 50 ' 51 In line with this model, TMS delivered to 
the left temporoparietal cortex possibly reduces AVH. 6,35,52 
Recently, we published a case of a schizophrenia patient with 
medication-resistant AVH that was successfully stimulated at 
the same region with transcranial direct current stimulation 4 
These findings suggest that AVH require the activation of 
speech perception neurocircuitry, consistent with what has 
been proposed as a model of underconstrained perception. 53 
Interestingly, none of the covariates demonstrated any 
prediction power, suggesting that neurobiological factors are 
more closely related to the possibility of modulating AVH with 
TMS. In addition, there was no statistical difference between 



the efficacy of theta-burst and 1 -Hz TMS protocols, which is in 
line with our recent finding that theta burst and 1-Hz TMS are 
equally effective in the treatment of AVH 24 

Several limitations of our study design merit comment. One 
must consider that no direct relationship between brain 
activity and the occurrence of hallucinations during scanning 
was investigated, because in ASL the mean cerebral 
perfusion over the scanning time is calculated and therefore 
temporal resolution is low. In addition, these results will have 
to be replicated in a sham-controlled design that applies the 
cutoff point found in our study. Both sham and stimulation 
groups should include an even number of patients with a high- 
enough level of STG perfusion before treatment. Thus, a 
decision-making algorithm for future studies could be devel- 
oped following the ROC analysis. A level of normalized CBF in 
the left STG of at least 0.26 could then be used to decide 
whether to include or exclude a patient because of expected 
treatment response in TMS. 

How can the higher STG perfusion in responders be 
explained? Given the fact that the target region for TMS was 
also situated in the STG, one might speculate that the efficacy 
of TMS relies on a certain level of brain activity in the area 
close to the target region. 29 This is in line with the suggestion 
that the disturbance of the AVH-generating network by TMS is 
the key factor 24 In the condition of relatively low brain activity 
in patients, an inhibitive TMS treatment might not be the 
appropriate treatment approach. Future studies should be 
aimed at investigating whether there are alternative brain 
regions in such nonresponders that could be involved in the 
generation of AVH, and that could be considered the possible 
target regions for TMS. 

In conclusion, this study showed that only patients with 
higher STG perfusion benefited from TMS as adjuvant 
treatment. Being noninvasive, relatively easy to use and 
independent of the occurrence of AVH during scanning, the 
measurement of resting perfusion with ASL before treatment 
might be a clinically relevant way to identify possible 
responders and nonresponders to TMS. 
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